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Concatenated Composite Pulses (CCCPs) are derived from various composite pulses widely employed in NMR and have been 
developed as high-precision unitary operations in Quantum Information Processing (QIP). CCCPs are robust against two systematic 
errors, pulse-length and off-resonance errors, in NMR simultaneously. We show experiments that demonstrate CCCPs are powerful 
and versatile tools not only in QIP but also in NMR measurements. 
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1. Introduction 

Nuclear Magnetic Resonance (NMR) is widely employed 
from Magnetic Resonance Imaging (MRI) in hospitals to pre¬ 
cision chemical analysis of various molecules in pharmaceu¬ 
tical companies |[Tj] thanks to highly developed NMR tech¬ 
niques. Some of these advanced techniques in NMR have been 
transferred to Quantum Information Processing (QIP) fH] since 
NMR manipulations are regarded as controlling and measuring 
quantum objects, called spins. 

We have been working on transferringone of existing NMR 
techniques, a composite pulse II II SB, to QIP and devel¬ 
oped high-precision unitary operations 171]. These operations 
are called Concatenated Composite Pulses (CCCPs) and are ro¬ 
bust against two systematic errors simultaneously at the cost of 
operation time, namely, the total length of pulses. These errors 
correspond to pulse-length errors (PLE) and off-resonance er¬ 
rors (ORE) in NMR manipulations. Instead of composing such 
pulses from scratch, CCCPs are designed by concatenating two 
composite pulses with different characteristics against system¬ 
atic errors. Thus these are descendants of existing composite 
pulses widely employed in NMR. QIP requires high-precision 
control of spins. Therefore, once CCCPs developed in QIP are 
employed for NMR measurements, significant improvement of 
signal strength is expected without any changes in the hardware 
settings. 
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The purpose of this paper is to feedback our achievement ob¬ 
tained in QIP to NMR measurements. In Sec. [2] we review the 
idea and properties of CCCPs. Sectionals the main part of this 
paper, where the superiority of CCCPs over the conventional 
methods in NMR measurements, single pulses and COSY, are 
demonstrated experimentally. Section[4]is devoted to summary 
and discussions. Some theoretical details are described with 
simulations in Appendix. 

2. Concatenated composite pulse (CCCP) 

In this section, we summarize the principle of CCCPs, i.e., 
composite pulses simultaneously robust against two typical sys¬ 
tematic errors, PLE and ORE, in NMR B5D . 

2.1. Pulse length errors (PLEs) and off-resonance errors 

(OREs) 

The system we consider throughout this paper is a nucleus 
with spin 1/2 (in short, a spin ) in a static magnetic field along 
the z-axis. We consider a rotation operation of the spin 

R(8, tp) = exp[-i 8n(cp) ■ cr/2\ (1) 

without errors. 8 is the rotation angle, «(</>) = (cos </>, sin </>, 0) is 
the rotation axis in the xy-plane, and cr = (cr x , cr y , erf) are the 
Pauli matrices. This rotation is generated by an ideal rf pulse in 
NMR. 

We consider a realistic square pulse in which both PLE and 
ORE are present simultaneously. We ignore the second and 
higher order error terms for simplicity. The rotation associated 
with a square pulse under only PLE is given as 

R' e ( 8, tp) = exp[-i(l + e)8n ■ cr/2] 

= R{8, (p) - i s8(n ■ (t)R{ 8, <p)t 2, (2) 
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where e is the strength of PLE, which is unknown but we as¬ 
sume that it is constant and small. Higer order terms beyond 
first order in e are suppressed in the second equality. This type 
of errors often appear because of inhomogeneity in the B\ field 
or lack of long term stability of spectrometers. On the other 
hand, the rotation associated with a square pulse under only 
ORE is given as 

R'f(6 , (/>) = exp[-i0(H • cr + fcr z )/2\ 

= R(9,<p) - \f sm(9/2)cr z , (3) 

where / is the strength of ORE. ORE is caused whenever the 
resonance of interest is not in resonance with the carrier fre¬ 
quency (transmitter frequency). In other words, ORE appears 
whenever there is a frequency offset between the carrier fre¬ 
quency and the Larmor frequency. Similarly to PLE, / is un¬ 
known but we assumed that it is constant and small. Therefore, 
when both PLE and ORE are present, the rotation associated 
with a square pulse is given as 

= exp[—i(l + e)9{n • tr + fcr z )/2\ 

= R{9,([>) - i s9(n ■ cr)R(0, </>)/2 - i/sin(0/2)crv. (4) 


BB1 



Figure 1: Schematic diagram of CCCP by using BB1 and CORPSE. Time goes 
from left to right in this pulse diagrams. 


as that in Eq. (0. We call the property residual error preserving 
property (REPP) with respect to ORE. 

In contrast to BB1, CORPSE is an IV = 3 composite pulse 
robust against ORE. CORPSE behaves as 

^corpse^’0) = R(9,<P) - i«(« ■ (r)R(9,4>)/2. (7) 


CORPSE possesses REPP with respect to PLE. 

We pointed out that not all the composite pulses have 
REPP |[I1. 


2.2. Composite pulses 

The NMR community has been developing a technique to 
overcome PLE and ORE by combining several pulses y), |5|, |6[]. 
Given a target rotation R(9 , 4>), we may find an equivalent rota¬ 
tion sequence which is equal to the target R(9, <f >) in case of no 
errors, as follows. 


R(9n, <Pn)R(9n-\ , <Pn-\) ' • -R(9i,<pi) - R(9,(f>). (5) 

Here, R(9,. (/>,) is the z-th rotation associated with the z-th square 
pulse and N denotes the number of square pulses. The point of 
the decomposition 0 is 

R'(9n,<Pn)R'(9n-i,<I>n-i) • • • R'(9\,(p \) + R’(9,(p ) 


if PLE and/or ORE exist. This is due to the non-commutativity 
among R(9j, </>,). Therefore, by appropriate tuning of the param¬ 
eters {9j,<f>j}^ = j in Eq. 0, we can obtain the sequence which i) 
virtually works as the target R{9, <p) when there are no errors, 
and ii) is less sensitive to the systematic errors. Indeed, such 
pulse sequences have been designed fisssaosn] by 
choosing the parameters [9i, <pj}f =1 in such a way that Eq. © 
has no first order terms with respect to the error parameter s or 
/. Note that these are the type A composite pulses, for which 
the first order cancellation applies independently of the initial 
state 

Now we present two typical composite pulses that are robust 
against either PLE or ORE: one is the Broad Band 1 (BB1) @] 
and the other is the Compensation for Off-Resonance with a 
Pulse SEquence (CORPSE) [jot], See, also Appendix A BB1 
under both PLE and ORE behaves as 


R’ ml (e,cf>) = R(9, </>) - if sin(9/2)<r z . (6) 

BB1 is obviously robust against PLE, but not against ORE. 
Moreover, the first order error term of BB1 is exactly the same 


2.3. Concatenated composite pulses (CCCPs) 

We show how to design a CCCP compensating both PLE and 
ORE simultaneously by taking the adva ntag e of REPP with 
BB1 and CORPSE, as an example [12, oT . BB1 is robust 
against PLE, while CORPSE is robust against ORE and has 
REPP with respect to PLE. Therefore, we replace all square 
pulses in BB1 with CORPSE (see Fig.0 This CCCP is called 
CORPSE in BB1, or CinBB in short. The number of square 
pulses in CinBB is 4 x 3 = 12. The number of square pulses 
in CinBB can be further reduced to N — 6 and the resulting 
CCCP is called the reduced CinBB, or R-CinBB for short. See 


Appendix A and [13] for further details. 


Other interesting approach to tackle both PLE and ORE was 
discussed by Jones 114]. He designed composite pulses com¬ 
pensating for higher order error terms of both PLE and ORE 
simultaneously. The rotation angle 9 is, however, fixed to n in 
these composite pulses. For the recent progress of the com¬ 
posite pulses including CCCPs, see the review by Merrill and 
Brown 11511 . 


3. Advantages of CCCP 

The advantages of CCCPs in NMR experiments are demon¬ 
strated with various experiments. We employ an ECA-500 
NMR spectrometer (JEOL RESONANCE Inc.). We evaluate 
the performance of the R-CinBB composite pulse as an exam¬ 
ple and compare the result with that of a square pulse. 

3.1. Single Pulses 

The experiments are carried out using 13 C-labeled chloro¬ 
form (Cambridge Isotopes) diluted in t/f,-acetone of approxi¬ 
mately 300 mM concentration at room temperature. In order to 
shorten T i for a rapid repetition, a relaxation agent (4 mM of 
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Figure 2: Spectra of FID signals after the 7r/2-square pulse (top) and the 
R-CinBB pulse (bottom) are applied as functions of PLE (s in Eq. 0). The 
strength of ORE is fixed to / ~ 0. 


Figure 3: Spectra of FID signals after the ^/2-square pulse (top) and the 
R-CinBB pulse (bottom) are applied as functions of ORE (/ in Eq.0. The 
strength of PLE is fixed to s ~ 0. 


Iron(III) acetylacetonate) is added. TVs and 7Vs are ~ 6 s and 
200 ms for 13 C, respectively, while those of proton are both ~ 
200 ms. We examine the performance of a composite 7r/2-pulse 
applied on 13 C and compare the result with that of a square 
pulse while proton is decoupled with WALTZ16 |Q. 

The determination of the pulse strength is not simple be¬ 
cause of inhomogeneity of B\ in the sample. We take a — of 
the 2;r-pulse duration as the nominal pulse duration of a square 
— pulse. It is obvious that the R-CinBB composite pulses are 
more advantageous compared with square pulses in terms of 
PLE, as shown in Fig[2] See, also numerical calculations shown 
in | Appendix B| and 11311 . 

We also examine the R-CinBB composite pulse in terms 
of ORE, as shown in Fig [3 The R-CinBB composite pulse 
is obviously more advantageous than the square pulse when 
-0.3 < / < 0.8 in Eq. Q. On the other hand, the spectra 
of the R-CinBB composite pulses corresponding to / < -0.5 
and 1.0 < / are more distorted than those of the square pulses. 
See, also numerical calculations shown in 
reference 0. 

3.2. COSY 

The two dimensional shift-COrrelation SpectroscopY 
(COSY) is one of the most important 2D-NMR measure¬ 


Appendix B and the 


ments Of We performe COSY experiments with 3-chloro- 
2,4,5,6-tetrafluoro-benzotrifluoride diluted in f/f,-benzen of an 
about 300 mM concentration. We observe signals of 19 F nuclei 
directly attached to a benzen ring ( 19 Fs at 2, 4, 5, and 6). TVs 
are between 0.6 and 1.0 s, while 7Ys are ~ 0.3 s. We use this 
molecule because of the following reasons. Firstly, as shown in 
Fig. SI 19 F signals of the molecule are widely spread. Secondly, 
the spectrum pattern is not so simple, although its molecular 
strcture is relatively simple. 

The pulse duration of a 7r/2-square pulse is 12.4 /rs which cor¬ 
responds to a Bj strength of 20 kHz in frequency unit. Note 
that the frequency difference between -160 and -118 ppm is 
20 kHz at 470 MHz and is comparable with the strength of B\. 
The total duration time of the R-CinBB pulse is 12.4 x 8 = 
99.2 /rs that is almost instantaneous like a square pulse in our 
experimental time scale and thus the replacement of a square 
pulse by the R-CinBB one leads no problem. The data are an¬ 
alyzed as follows. For both t\ and 12 period, shifted sine-bell 
window function is multiplied. For t \, zero-filling is done once. 
This time-domain data is then Fourier-transformed and phased. 

B 1 strength of a 7r/2-square pulse in frequency is compara¬ 
ble with the frequency difference between the highest and low¬ 
est peaks. Therefore, we expect that COSY experiments with 
square pulses may have some difficulties. There are attempts 
with BURP (Band-selective, Uniform Response, Pure Phase) 
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Figure 4: 19 F spectrum of 3-chloro-2,4,5,6-tetrafluoro-benzotrifluoride ride in 
d6-benzen. Each peak is enlarged to show detailed structures (1 ppm). Peak 
assignments are as follows; 1 (-118.0 ppm), 2 (-126.5 ppm), 3 (-135.5 ppm), 
and 4 (-160.5 ppm) are identified as F2, F4, F6, and F5, respectively. 19 F signal 
of trifluoromethyl group is not observed in this region. 


to overcome these difficulties 00. 

As we expected, the correlation peak between -118.0 ppm 
(fl) and -126.5 ppm (f2), and that between -126.5 ppm (fl) 
and -118.0 ppm (f2) are hardly visible in the case of square 
pulses. On the other hand, these are much higher intensities in 
the case of R-CinBB pulses, although they are small as is found 
in Fig. 0 The advantage of R-CinBB pulses are more clearly 
demonstrated in 1-D spectra, as shown in Fig. [6] The phases of 
peaks obtained with square pulses are highly distorted and this 
may be the reason why the above correlation peaks are hardly 
visible. 
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Figure 5: 2D-spectra obtained with (a) a square and (b) an R-CinBB pulses. 


4. Summary 

Composite pulses have been developed in the NMR commu¬ 
nity and are widely employed in daily measurements. Concate¬ 
nated composite pulses (CCCPs) are directly descended from 
them and have been developed as robust unitary operations for 
quantum information processing. They are robust against two 
systematic errors, the pulse length error and off-resonance error 
in NMR, at the cost of execution times. Note, however, that this 
cost is not important in usual experimental conditions of liquid- 
state NMR and thus CCCPs are always advantageous compared 
with square pulses. 
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Appendix A. Example of composite pulses 

We present details of the Broad Band 1 (BB1) j@], the 
Compensation for Off-Resonance with a Pulse SEquence 
(CORPSE) H and the reduced CORPSE in BB1 (R- 
CinBB) foil . 

Appendix A.l. BB1 

BB1 is an N — 4 composite pulse robust against PLE. The 
parameters are chosen as 

6\ — 63 — 1 r , O2 — 27 r , O4 = $, 

<Pi - 03 = f + arccos[— 0 /( 47 t)], <p 2 = 2>(f>\ - 2<p, $>4 = <p. 

(A.l) 

BB 1 under both PLE and ORE behaves as 

R'bbX 6 ' 0 ) = R '(0, <P) R 'fa <t>i)R'(2x, <p 2 )R'(n, cp 1 ) 

= R(fi, 0) - if sin(0/2)cr z . (A.2) 
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For instance, the sequence implements zr-pulse when 

01 = 0 3 = 05 = 7 r, 0 2 = 27 t, 

04 = 77t/3, 05 = 5tt/3, 06 = n/3, 

0i = 03 = tt - arccos(l/4), 02 = 30i, 

04 = 06 = o. (A.6) 


Appendix B. Simulations of NMR experiments 

We perform some numerical simulations in order to support 
experiments. 



5(19F) / ppm 

Figure 6: lD-spectra obtained with (a) square and (b) R-CinBB pulses in COSY 
measurements. 

Appendix A.2. CORPSE 

CORPSE is an N — 3 composite pulse robust against ORE. 
Its parameters are 

0i = 2«i7r + 0/2 — k, 02 — 2ri27i — 2k, 03 = 2nj,n + 0/2 - k, 

01 = 02 - 7T = 03 = 0, £ = arcsin[sin(0/2)/2], (A.3) 


Appendix B.l. method 

We have to take into account of a non-unitary time develop¬ 
ment caused by a spin-spin relaxation with a characteristic time 
72, in order to simulate NMR experiments. We introduce this 
effect as a phase flip channel Cl¬ 
in the case of single spin experiments, 

p(t + A) = p,JA)p(t) + (1 - p ss ( A)) Ad (cr z ,p(t)) , 

t w l+exp(-A /T 2 ) 

Pssi A) = ---, 

where Ad(*,p) = *p*'' and A is a small time interval. 

During a pulse, the time evolution is simulated as follows. 

p\t + T p ) = p ss (r p )p(i) + (l -pssiTp)) Ad (cr z ,p(t)), 
p(t + T p ) = Ad(U puhe ,p'(t + t p )), 

where U pu i se is a unitary operation generated by the pulse. Note 
that t p is the total pulse duration and is assumed to be small. 
Therefore, we employ the Suzuki-Trotter formula [20], here. 

In the case of COSY experiments, the simulations during the 
evolution and detection periods Cl are done as follows. 


where n\,n 2 ,n 3 are non-negative integers. In particular, when 
we take «i - nj =0 and n .2 - l, the execution time is mini¬ 
mized. Then CORPSE is called short CORPSE. Other notable 
case takes place when n\ - «2 + «3 = 0. In this case with both 
PLE and ORE, CORPSE behaves as 

^corpse^’ 0) = </>) R '( 0 2 ,0 + n ) R \ 6 \, 0) 

= R(9, 0) - i e(n ■ <r)R(9, 0)/2. (A.4) 

Appendix A. 3. reduced CORPSE in BB1 
The R-CinBB is given as follows. 

01 = 03 - tt, 02 = 2i r, 

04 — 06 + 2u — 2u + 0/2 — k, 05 = 2n — 2k, 

0i = 03 = 0 + arccos[- 0 /( 47 r)], 0 2 = 30i - 20, 

04 = 05 — n = 06 = 0, k = arcsin[sin(0/2)/2]. (A.5) 


p'(t + 6) - 


+ 


+ 


pit + 8) = 




Pit ) 


6 

—— Ad(cr- ® cr o ,p(0) 

Zi 2,1 

5 

—— Ad(o- 0 ®cr,,p(f)), 

^7 2,2 

cr- ® cr, 

Ad(exp(-74^—j |—-),p it + 5)), 


where 72,; is the T 2 of the 7-th spin, pit + nS) can be obtained 
by iterating the above operations n-times. We assume that the 
scalar coupling is in a weak coupling limit 12ill . Note that 
the above decomposition is based on the Suzuki-Trotter for¬ 
mula 1 20!l by assuming 5 is sufficiently small and that p ss (6) ~ 
1 - d/(274 ,j for the 7-th spin. 

During a pulse, the time development is similarly simulated 
as in the case of single spin experiments. 


5 





























Appendix B.2. Single pulses 


The signal amplitude is estimated in the case of a typical 
NMR measurement where a control field (B\ in NMR) is inho¬ 
mogeneous over the sample and assumed to have the Gaussian 
distribution |22|. 


P(Bi) 


1 

y[2ncr 


exp 


CB 1 -.B 1) 2 

2<x 2 


(B.l) 


where B\ is the average strength of B\ and cr determines the 
inhomogeneity. On this assumption with a previously esti¬ 
mated value of cr = 0.05 1 221 ] with our NMR spectrometer, 
the signal strength distributions after a nominal 7r/2-pulse (or, 
B\t p = 7t/2 where t p is the pulse duration) is calculated, as 
shown in Fig. IbTTI It is obvious that R-CinBB is advantageous 
over the square pulse. 



position 


Figure B.7: (colour online) Projection of magnetization on the xy-plane after a 
pi/2-pulse. The red (blue) curve is the outcome of a square pulse (the R-CinBB 
pulse) as a function of z along the sample. The gray line indicates the strength 
of B] as a function of z. 


A phase error in the NMR spectrum occurs under ORE, as 
shown in Fig. IB. 81 Here, / = 0.1 is assumed. The red curve 
in this figure is a simulated result of a square 7r/2-pulse and ap¬ 
pears distorted. On the other hand, the blue curve is a simulated 
result of a R-CinBB-w/2-pulse and is not distorted. 



20 30 40 50 60 70 80 90 100 110 


frequency 


Figure B.8: (colour online) Simulated NMR spectra under ORE. The red curve 
is a simulated result of the square /r/2-pu]se with / = 0.1 and appears distorted. 
On the other hand, the blue curve is a simulated result of R-CinBB pulse with 
the same / and is not distorted. 


Appendix B.3. Hahn echo 


Two pulses in Hahn echo experiments [2j] are affected by 
PLE and ORE, thus measured 72 under such errors may be er¬ 
roneous. We simulate an NMR equipment with fluctuating PLE 
and ORE by varying PLE and ORE in every spin echo measure¬ 
ment, as shown in Fig. IB.91 Their means are s = f = 0.1 and 
their standard deviations are both 0.08. Although these figures 
may be unreasonably large for modern NMR spectrometers, we 
can obtain the correct 72 with R-CinBB pulses. 



Figure B.9: (colour online) Semi-log plot of echo signal with square pulses (red 
line) and R-CinBB pulses (blue line) as a function of the waiting time r. Black 
dashed line is in error-free case. 


Numerical simulations of the Hahn echo experiments as a 
function of the error strengths are summarized in Fig. IB. 101 
The Hahn echo experiments with square pulses (left figure) are 
strongly affected by PLE, while those with R-CinBB pulses are 
robust against these errors. It turns out that a composite pulse 
robust against PLE is sufficient in order to obtain a correct 72 
under both PLE and ORE. 



Figure B.10: (colour online) Simulated 72 with square pulses (left) and with 
R-CinBB pulses as a function of PLE and ORE. The Hahn echo experiments 
with R-CinBB pulses give us correct 7Vs even in erroneous cases. 


Appendix B.4. COSY 

Here, we simulate COSY experiments applied for two-spin 
molecules of which interaction is a scalar coupling in a weak 
coupling limit, as shown in Fig IB. 1 ll The chemical shifts of 
these spins are assumed to be 1 and 4 ppm and J = 0.5 ppm. 
In COSY experiments, spurious peaks called axial peaks some¬ 
times appear due to inaccuracy of the first pulse yj] . We are 
able to reproduce these axial peaks in the simulation of COSY 
experiments with square pulses (e = / = 0.1), as shown in 
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Fig. IB. 1 ll On the other hand, no axial peaks appear in the sim¬ 
ulation with R-CinBB pulses. 


Original COSY 


CCCP-modified COSY (R-CinBB) 



Figure B.ll: (colour online) Simulation of COSY experiments of a two-spin 
molecule with with square pulses (left) and R-CinBB pulses (right), where e = 
/ = 0.1, J = 0.5 ppm and chemical shifts are 1 and 4 ppm. Spurious axial 
peaks (green peaks) are observed in the simulation with square pulses, while no 
such peaks are observed in that with R-CinBB pulses. 
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